Wave-packet Formalism of Full Counting Statistics 
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We make use of the first-quantized wave-packet formulation of the full counting statistics to 
describe charge transport of noninteracting electrons in a mesoscopic device. We derive various 
expressions for the characteristic function generating the full counting statistics, accounting for 
both energy and time dependence in the scattering process and including exchange effects due to 
finite overlap of the incoming wave packets. We apply our results to describe the generic statistical 
properties of a two-fermion scattering event and find, among other features, sub-binomial statistics 
for nonentangled incoming states (Slater rank 1), while entangled states (Slater rank 2) may generate 
super-binomial (and even super-Poissonian) noise, a feature that can be used as a spin singlet-triplet 
detector. Another application is concerned with the constant-voltage case, where we generalize the 
original result of Levitov-Lesovik to account for energy-dependent scattering and finite measurement 
time, including short time measurements, where Pauli blocking becomes important. 

PACS numbers: 73.23.-b, 73.63.Nm, 73.50.Bk, 05.60.Gg 



I. INTRODUCTION 

Charge transport across an obstacle in a wire is a 
statistical process, whose complete description is pro- 
vided by the probability function P(n,t), telling how 
many charge carriers n are transmitted through the wire 
during the time t. The calculation of this full count- 
ing statistics usually aims at the generating function 
X(A, t) = P(n, t) e lXn for this process, from which the 
probability distribution P(n, t) follows through simple 
Fourier transformation ,F[x(A, t)] = P(n,t). The proper 
physical definition of the generating function x(A, t) is a 
nontrivial problem and has been solved by Levitov and 
Lesovik back in 1993,— see also Ref. [H, with numerous 
applications to follow^. The original definition includes 
a 'charge counter' in the form of a spin, coupled via the 
gauge potential to the moving charges, and has been cast 
in a second-quantized formalism of appreciable complex- 
ity. The recent observation^ of the correspondence be- 
tween the generating function Xi(A) of the full counting 
statistics for one particle and the notion of fidelity in a 
(one-particle, chaotic) quantum system^ has lead to a 
much simpler first-quantized formulation of full counting 
statistics, including the generalization Xiv(A) to N parti- 
cles. In fact, a first-quantized version of charge transport 
to calculate noise has been already introduced some years 
ago£ Furthermore, such a wave-packet formalism natu- 
rally describes the statistics of pulsed transport, where 
unit-flux voltage pulses generate single-particle excita- 
tions feeding the device of interesl! 2 -^^^ (a source in- 
jecting individual electrons into a quantum wire has been 
realized in a recent experiment-^ 1 ) . The simplicity of the 
first-quantized formalism then has allowed to obtain non- 
trivial results on the full counting statistics for an energy 
dependent scatterer, including its dependence on the ex- 
change symmetry of the transported charge.— 

In this paper, we make intense use of this wave-packet 



formalism of charge transport and (re-)derive various ex- 
pressions for the characteristic function x./v(A) in a much 
simplified manner. We start with an TV-particle Slater 
determinant made from orthonormalized single-particle 
wave functions <p m describing fermions incident from the 
left and derive the associated characteristic function de- 
scribing the full counting statistics in determinant form, 



Xiv(A) = det(^ m |l -T + Te lX \<j> n ) 



(1) 



with the operator T describing the energy dependent 
transmission across the scatterer, T — J (dk / 2n)Tk\k) (k\ 
in momentum (fc) representation (here, the particle 
number N replaces the time variable t in the original 
formula^). The determinant in Eq. fl} can be cast in a 
product form 
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Xiv(A)= 



(2) 



where r m are the eigenvalues of the Hermitian operator T 
in the space spanned by the basis states | <j> n ). We denote 
the distribution in ([2|) as 'generalized binomial'. 

In a real experiment, the unit-flux voltage pulses gen- 
crating the incoming wave packets may overlap. For this 
situation, we rederive the simple and elegant expression 
([2|) for the full counting statistics, but with the coef- 
ficients r m now replaced by the roots of a generalized 
eigenvalue problem incorporating all effects of fermionic 
statistics and the full energy dependence of the trans- 
mission. The results (JI]) and ([2|) apply to a nonentangled 
incident state in the form of a Slater determinant^ an 
extension to include entangled states of Slater rank 2 is 
provided as wellp^ 

Next, we generalize the result (JlJ to describe a setup 
where both the scattering process and the counting win- 
dow depend on time and find a compact result in the 
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form ((T|) with 

T^T Q =U^QU, (3) 

where IA denotes the single-particle time evolution oper- 
ator and the operator Q projects the wave function onto 
its measured (counted) part. Full counting statistics for 
fcrmionic atoms in determinant form has been derived in 
Ref. UM through transcription of the bosonic expression^ 
to the fermionic case, see Ref. Il7l for a recent application. 

Finally, we extend the result Eq. ^ to the situation 
where the incoming state consists of an incoherent su- 
perposition of many Slater determinants with different 
particle numbers. For the case of particles incident only 
from the left side, we find the result |T]) with 

T^rfT Q , (4) 

where rj denotes the one-particle occupation-number op- 
erator. In addition, the determinant in (TT]) has to be 
taken over all the single-particle Hilbert space. 

We make extended use of these formulas: for a two- 
particle problem, we show that, i) an incoming state 
described via a simple Slater determinant cannot gen- 
erate a Fano factor F = ((n 2 )} / (n) > 1 — (n)/2 (i.e., 
noise is always sub-binomial and in particular also sub- 
Poissonian; there is no bunching); the above cumulants 
are obtained through the generating function x(A) ym 
= (-i) d{logx\\=o, n) upon proper choice of T k , 
an entangled incoming state can generate any value for 
the Fano factor F < 2, and Hi) for two spin- 1/2 fermions, 
we show that a simple scattering experiment provides in- 
formation on the entanglement of the incoming state (cf. 
also Ref. [II). 

Subsequently, we analyze the situation with N 
fermions and derive the full counting statistics for a con- 
stant voltage (V) drive, thereby generalizing the original 
result of Levitov and Lesovik^ to describe transport with 
an energy dependent scattering transmission (cf. Ref. [l9l) . 
Our result, 

logx^A) = N 2 -^ / J* log(l -T k + lie"), 
eV J 2tt 

admits the simple interpretation of the full counting 
statistics as deriving from the transmission of the un- 
balanced Fermi sea residing between energies E F and 
E F + eV, with E F denoting the Fermi energy and V the 
applied bias. Using an alternative derivation based on 
([3|) and stationary scattering states, we determine the 
short-time limit of the counting statistics and rederive 
the binomial result (0 in the long-time situation, with 
the particle number N replaced by the measuring time t, 
N — ► t eV/2irh. The use of our determinant formula com- 
bined with Szego's theore m 20 i 21 will allow us to present 
a rigorous derivation of these results. 

In the following, we give a short review of previous 
work on the subject and then derive the characteristic 
functions (HJ and (0 of N incoming fermions. In Sec. lIIIi 



we apply these results to discuss the statistical transport- 
properties of two fermions. Section ITVl is devoted to the 
calculation of the characteristic function for the constant 
voltage case starting from iV-particle trains and letting 
the width of the individual wave packets go to infinity. 
In Sec. El we derive the results ([3J and (0J describing the 
setup involving a time dependent scattering and count- 
ing incoherent superpositions of incoming particles. We 
rederive the constant voltage result as an application, in- 
cluding the short-time limit. 



II. FULL COUNTING STATISTICS 

The first suggestion^ 2 , of a generating function for full 
counting statistics relied on the straightforward expres- 
sion x(A,i) = (exp[iA J dt'X(t')]), wherel(i) denotes the 
current operator. It then was soon realized^ that this 
definition does not correspond to any known (even on 
the level of a 'Gedanken Experiment') measuring pro- 
cedure; still, this first definition produced the correct 
results for all irreducible zero-frequency current-current 
correlators ({Xq ■ ■ - To)}) (see also the discussion in Ref. 
[23h . The first 'practical' definitional of a generating func- 
tion x(A, t), corresponding (at least in principle) to a real- 
istic counting experiment, involved a spin-galvanometer 
as a measurement device (see also Ref. y). Recently, 
it has been pointed out! that this suggestion (corre- 
sponding rather to a 'Gedanken Experiment') could ac- 
tually be realized with qubits serving as a measuring de- 
vice, whereby the 'environmental noise' generated by the 
transmitted charge serves as the measurement signal for 
the full counting statistics. This contrasts with the usual 
interpretation of the 'environmental noise' as being re- 
sponsible for the qubit's dephasing 2 ^ expressed through 
the fidelity, and also relates to the competition between 
the gain of information and dephasing— in quantum mea- 
surement theory. 

The insight on the equivalence between the notions of 
fidelity and full counting statistics has motivated a first- 
quantized formalism of the counting problem in terms of 
wave packets. Fidelity |xfid|, the modulus of the overlap 
Xfid = (^2 1*1)1 was introduces by Peres^ in the context 
of chaotic systems. It measures the overlap between two 
wave functions 1^1,2) which describe an initial state |^o) 
which has evolved under the action of two slightly differ- 
ent Hamiltonians. In the context of full counting statis- 
tics of a single particle measured by a spin counter, the 
wave functions \&i and "J 2 are substituted by scattering 
states ^out an d 'I' out interacting with the spin counter 
in the states |f) and ||), resulting in an expression for 
the generating function in the form \i — (^outl^out)- 
This new first-quantized formulation in terms of wave- 
packets provides a drastic simplification as compared to 
the original second-quantized formalism. 2 While the use 
of a second-quantized formalism is mandatory for the 
description of particles describing bosonic excitations of 
fields (photons, phonons, etc.), here, we deal with non- 
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relativistic electrons where the particle number is fixed, 
thus allowing for an alternative first-quantized descrip- 
tion. Moreover, our wave-packet formalism has technical 
merits (e.g., in the description of energy dependent scat- 
tering or in the classification of two-particle scattering 
events) and also provides a better physical understand- 
ing. We remark, however, that in dealing with finite tem- 
peratures we make use of the second-quantized formalism 
in Fock space. 

An alternative method, to the procedure based on a 
spin counter, was pursued in several contributions 27,28 - 29 
where the full counting statistics and, in particular, its 
generating function x(A, i), was constructed using only 
basic quantum mechanical definitions; starting with an 
initial state in the form of an eigenstate of the particle 
number operator with a fixed particle number to the right 
of the scatterer (or the 'counter'), a second projection (to 
eigenstates of the number operator) onto the final state 
is carried out after the observation time t. Both pro- 
cedures, projection and spin-counting, lead to the same 
expressions for the generating function x, provided that 
the incoming state involves no superposition across the 
scatterer. In the latter situation, the explicit calculation 
using a spin-counter produces a fidelity describing the de- 
coherence of the spin, while an interpretation in terms of 
a generating function can produce probabilities for non- 
integer charge transport^ and hence is unphysical. On 
the other hand, the projection method, destroying such a 
superposition in the course of the first measurement, al- 
ways admits an interpretation in terms of probabilities. 



A. One particle 

In this paper, we make extensive use of the first- 
quantized formulation of the generating function: start- 
ing with a simple one-particle problem, we exploit 
the equivalence between the notion of fidelity and full 
counting statistics^ Consider an incoming wave packet 
ip(x; t — » — oo) from the left of the form 



dk 
2^ 



4>i(k)e 



ikx — ie(k)t 



(6) 



with normalization J(dk/2n)\(j)i(k)\' 2 = 1, cf. Fig. [TJ In 
the following, we assume (for simplicity) a linear spec- 
trum e = v F k with v F the Fermi velocity; at low temper- 
atures and voltages the interesting physics usually takes 
place near the Fermi points. The momentum Hk and the 
energy he are measured with respect to the Fermi mo- 
mentum k F and the Fermi energy E F . Here and below, 
the wave-packets include only momenta with k > in 
order not to disturb the Fermi sea which is considered to 
be the vacuum in our analysis. The scatterer at x — 
is characterized by momentum(energy)-dependent trans- 
mission (reflection) amplitudes t k (r^; particle reflection 
takes us to the branch e = —v F k, with k measured rela- 
tive to —k F ). The spin- (or qubit-) counter, placed to the 
right of the scatterer, contributes a phase-factor e ±lA / 2 



to the wave function, where the sign depends on the state 
| J.) of the spin-counter. The outgoing (t — + oo) wave- 
function assumes the form (we place the counter right 
behind the scatterer at x = 0) 



dk 
2^ 



r k e- lk(x+VFt) e{-x) 



(7) 



and consists of reflected (x < 0) and transmitted (a: > 0) 
parts; Q(x) is the unit step-function. The fidelity Xi(A) 
is given by the overlap of wave functions with slightly 
different perturbations in their evolution, here, with cou- 
pling to opposite spin-configurations |f) and ||), 

= J dxip- ut (x;t)* iP+ ut {x;t) 

( "° 0) f ^(l-T k+ T k e^)\Mk)\ 2 



II -T + Te 



iA I 



(8) 



in the asymptotic or long-time limit, the integration over 
space is trivially done by exploiting the complete sep- 
aration of the wave function into transmitted and re- 
flected parts. Furthermore, the time dependence disap- 
pears as soon as the transmitted wave function has passed 
the counter. The transmission probabilities T k = \t k \ 2 
are the eigenvalues of the transmission operator T = 
J (dk/2n)Tk\k) (k\. Given the above specific coupling to a 
spin, the fidelity is equivalent to the characteristic func- 
tion 



X 



(A)=^P m e l 



(9) 



of the full counting statistics as defined in Ref. 24, where 
a spin-galvanometer has been used as a measuring de- 
vice. The Fourier-coefficients P m are the probabilities for 
transmitting m particles. For the simple example of one 
incoming particle only two outcomes are possible, parti- 
cle reflection with probability Po = 1 — (T) and particle 
transmission with Pi = (T), where (T) — (</>i|T|0i) de- 
notes the average transmission probability. Knowing the 
characteristic function, the cumulants ((n- 7 )) can be ob- 
tained as the coefficients in the Taylor series of logx(A), 



{{n3)) = (Ay ' iogx(A) 



A=0 



(10) 



The ratio F = ((n 2 ))/(n) between the second and the 
first cumulant, called Fano-factor, will be of special in- 
terest later. 



B. N particles 

Next, we extend the above description to N particles 
with an incoming wave function ^E'(k) defined in momen- 
tum space; the vector k = (fci, . . . , kisr) specifies the N 
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fi(k) / 2 (fc) 




Ej? J 


eV(t)6 



Xy 



X s = 



FIG. 1: Quantum wire with scattering center located at x a 
giving rise to a momentum-dependent scattering probability 
T k . A time-dependent potential eV(t) applied at xv (to the 
left of the scatterer) generates incoming wave packets /i, /2 
with overlap S = (/2I/1). A counter, placed at x c (to the 
right of the scatter) , measures the statistics of the number n 
of transmitted particles. In our analysis, we consider incom- 
ing wave packets with momenta k > residing outside the 
Fermi sea. As a result, the Fermi sea, which is not accounted 
for in our analysis, is not disturbed in the asymptotic time 
limit. For finite measuring times, the presence of the Fermi 
sea generates additional equilibrium noise which we do not 
consider in this article. 



momenta of the particles. We consider independent par- 
ticles without interaction which scatter independently. 
After scattering, the outgoing wave function assumes the 
asymptotic (t — ► 00) form 

l/&t(*5*) = { ft [r km e- ik ^ Mt M-Xm) 

m=l 

+ t km e* k ^ x ™- v ^e ±lX / 2 Q(x m )}}y(k), (11) 

i.e., the evolution is the product of the single- 
particle evolutions in expression ([7]). The character- 
istic function of the full counting statistics xjv(A) = 
J <ix -0~ ut (x; i)*'i/)^ ut (x; t) then can be cast into the form 



XN 



( A ) = {n/ ^(i-T fcm+ r feme ^)}|*(k)| 2 . (12) 



So far, we did not specify the specific type of incoming 
wave function. If we limit ourselves to Slater determinant 
states composed of orthonormalized single particle states 



*(fci,...,fcjv) = -~,= det0 m (fc n ), (13) 

the expression Eq. IT2|) can be rewritten as a single de- 
terminant (see Eq. (|rj5|) ) 

/dk 
— ^ m (k)(l-T k + T k e lX )Mk) 

= det(0 m |l-T + 2V A |0„) (14) 

involving the single-particle matrix elements ((f> m \0\(f) n ) 
of the operator O = 1 - T + Te lX . 



C. Nonorthogonal basis 

In a physical realization of such a scattering experi- 
ment, one usually does not populate orthogonal states 
as used in the above construction of the Slater determi- 
nant. E.g., in the setup of Fig. [T] the electrons typically 
occupy states /1 and /2 with a finite overlap, i.e., they 
are nonorthogonal. Of course, an A^-particle Slater deter- 
minant can be constructed as well out of nonorthogonal 
states \f m ), provided they are linearly independent, i.e., 
det(/ m |/„) 7^ 0. The properly antisymmetrized and nor- 
malized wave function (|13p then acquires the form 



1 



^iV!det(/ m |/„ 



:det/ m (fc„). (15) 



^ f (h,...,k N ) = 



Inserting this expression into (fl2j) and repeating the cal- 
culation that led to (fl~4|) . we obtain the generating func- 
tion in the form of a ratio of two determinants, 



Xiv(A) = 



det (/ m 1 1 - T + Te lX \ f n ) 

det(/ m |/„) 
det(S- f - V + T f e iX ) 



det Sf 

with the two N x N matrices 

S™„ = (fmlfn) and T£,„ = 



(fm\T\f n 



(16) 



(17) 



D. 



Invariance of Slater Determinants under Linear 
Transfer mat ions 



It turns out that the expression (T6")) for the generating 
function can be drastically simplified and rewritten in a 
generalized binomial form. As a first step towards this 
goal, one has to realize that an A^-dimcnsional Hilbcrt 
space Hn, spanned by the single-particle wave- functions 
fn(k) 7 defines exactly one properly antisymmetrized wave 
function, or, equivalently, there exists (up to a phase fac- 
tor) only one associated A^-particle Slater determinant 
state. The antisymmetrized A^-particle state is thus a 
property of the Hilbert space Hjy and is independent on 
the basis chosen^ 

Consider, as a simple example, a two-particle Slater- 
determinant state (in second-quantized notation) \^>) = 
a 2 a I|0)j with the vacuum-state |0) and Fermionic opera- 
tors ai,2- Defining the new operators a± — (ai±a2)/v / 2, 
we easily see that the two-particle state 

«l a -|0> - 5(4 + 4)(«I-4)|0) =4a\\0) = |*> (18) 

remains unchanged. Consider then a general iV-particle 
Slater determinant state of the form Eq. (fTS")) . Trans- 
forming the basis states f m (k) to new states g m (k) via 
the complex linear transformation 



9m{k) 



]Kmfn(k), detA^O, 



(19) 



5 



the antisymmetric combination 

detg m (k n ) = (detA) det/ m (fc n ) (20) 

remains invariant up to the factor det A; here, we have 
used the fact that the determinant of the product of 
two matrices is the product of the individual determi- 
nants. Furthermore, the normalized iV-particle Slater- 
determinant states ^ and \I/ 9 obey the relation 

...,k N ) = sgn(det A)V f (h, ...,k N ) (21) 

with sgn(a;) = ac/| as j . The only effect of adopting a 
new basis is the appearance of an overall phase factor 
sgn(det A) which drops out of the characteristic function 
(fT2"]) . Therefore, the full counting statistics calculated in 
the bases / and g give identical results. 

E. Diagonalization 

The above invariance can be used to simplify the cal- 
culation of the full counting statistics. Furthermore, even 
without specification of the (time- independent) scatterer, 
one can obtain valuable insights about the structure of 
possible outcomes in the counting statistics. In partic- 
ular, it turns out that the most general full counting 
statistics for a Slater-determinant state is given by a gen- 
eralized binomial expression of the form ([2]) . 

Let us first investigate how the invariance under linear 
transformations, Eq. (fT9|) . manifests itself in the deter- 
minant formula Eq. Q16p . To this end, we note that any 
single-particle matrix B of the form (| 1 T[) transforms un- 
der the linear transformation A of the basis functions 
according to 

B 9 =A t B / A, B = S,T. (22) 

Since det(AB) = det A det B, we find that the character- 
istic function \n (we define X-^ = — 4- T^e lA ) is 
invariant under the change of basis, 

_ dctX^' _ | det A | 2 det X^ _ detX 9 
XN ~ detS/ ~ | det A | 2 det Sf ~ detS^' ( ' 

This invariance can be exploited by going over to new 
orthogonal basis functions g m (k) with an overlap matrix 
S 9 m = S mn and a transmission matrix assuming a diago- 
nal form T 9 nn = T m 5 mn . The possibility of simultaneous 
diagonalization of the matrices T 9 ^ and is a con- 
sequence of the transformation law (f2"2")l characteristic of 
bilinear forms (as opposed to linear transformations L 
which transform according to L 9 = A _1 L^A) combined 
with the positivity of S* . The corresponding eigenbasis 
g m and eigenvalues r m of T 9 ^ can be found by solving 
the generalized eigenvalue problem 

{T f - T m S f )a rn = (24) 



with the normalization aJ rl S^a m = li^L The eigenvec- 
tors a m constitute the column vectors of the trans- 
formation matrix A = (oi, . . . , djv)- The eigenvalues 
are given by the roots of the characteristic polynomial 
det(T / -t S f ) = 0. The full counting statistics, Eq. (fl6|) . 
written in the new basis g m (k) assumes the generalized 
binomial form 

N 

Xiv(A) = ]J(l-r m +T m e iX ), (25) 

m— 1 

where the determinant has been evaluated explicitly and 
the result depends only on the eigenvalues r m . The gen- 
eralized eigenvalue problem can be reduced to a nor- 
mal one by rewriting the problem in a orthogonalized 
basis (j) m (k), with = Ijv, which can be obtained 
by the Gram-Schmidt procedure or by setting <f> m {k) = 

EJ(S / )- 1/2 ]« m /n(fc). 

From the above, we see that the concrete form of 
the eigenvalue problem (|2"4"|) is basis dependent, whereas 
the eigenvalues and vectors are simply a property of the 
transmission operator T operating in the Hilbcrt space 
Hn with the scalar product (f\g). Indeed, it is possible 
to find the eigenvalues and eigenvectors in a basis inde- 
pendent way using the positive definite quadratic forms 
T (g) = (9\T\g) and S(g) = (g\g), g S H N . Representing 
the bilinear form T{g) with fixed S(g) = 1 as a polar 
plot with T(g) the radius and g defining the direction 
in Hn , we obtain an ellipsoid in A^-dimensional space. 
The lengths of the main axes of this ellipsoid then con- 
stitute the eigenvalues and the associated directions the 
eigenvectors of the problem (pMjl i^ 2 - The eigenvalues r m 
are constrained to the interval [0, 1] as T(g) > and 
T{g) < S{g) due to unitarity. 



F. Full Counting Statistics for Entangled States 

The above discussion has concentrated on incoming 
states described by a single Slater determinant, i.e., 
nonentangled states with Slater rank 1. It is instructive 
to generalize this discussion to entangled states involv- 
ing a coherent superposition of Slater determinants. We 
start from an incoming state of N particles with Slater 
rank 2, 

*(k) = a* I (k) + /3^ II (k), (26) 

where \I' I (k) and ^ n (k) are normalized iV-particle Slater 
determinants describing particles incoming from the left 
and made from single particle states fm(k) and fm{k), 
m = 1, . . . , N\ the complex numbers a and /3 have been 
chosen such as to make 'l'(k) normalized. The character- 
istic function for the full counting statistics (TT2"]) assumes 
the form 



6 



XN 



N 

w=[n 



m— 1 



u,n m 

~2T 



(1 - T km + T fcm e lA ) laH^Ck)! 2 + |/3| 2 |* n (k)| 2 + 2Re{a/3** I (k)* II (k)*} 



(27) 



r 



where Re denotes the real part. The first two terms re- 
duce to generating functions for simple Slater determi- 
nant states and we can write 



X ,v(A) = |a| 2 X ] v (A) + |/3|V^(A) 

+ Q rxf(A)+a*^f(-A)" 



(28) 



with 



x5v(A) 
X"(A) 
X?T(A) 



det(S^ -T^ 1 +T^e iA ) 

detS/ 1 
d et(S/ n -T/ n +T/ n e ' ;A ) 

detS/ 11 
det(S mix - T mix + T mix e iA ) 

VdctS/'S/ 11 



(29) 



The matrices with superscripts f l and f u have been de- 
fined in Eq. (jTTJ) , while the new Hermitian matrices with 
a superscript 'mix' are given by the mixed matrix ele- 
ments 



cmix 

mn 



The first two terms in 
cf. 



(Ofi), T^ = (f^\T\f n ). (30) 
can be diagonalized as before, 



N 



Xiv(A)= Xlil-rl+Tle^), (31) 

m— 1 
N 



with the eigenvalues and given by the roots of 
det(T : - ^S 1 ) = and det(T n - r n S n ) = 0. 

Let us then concentrate on the characteristic function 
X mlx (A). Unfortunately, there is no generic procedure to 
follow in this case, as the matrices S mlx and T mlx are not 
Hermitian any more and hence the expression (|27p can- 
not be further simplified in general. In particular, the 
characteristic function x mlx (A) is not invariant under in- 
dividual transformations of the bases and (such 
basis transformations leave the Slater-determinants in- 
variant only up to a phase factor, which dropped out in 
the calculation of the characteristic function of a single 
Slater determinant state but does not when two Slater 
determinants are superimposed coherently). In order 
to proceed further, we restrict ourselves to specific sit- 
uations where S mix = or detS mix ^ 0. The most 
trivial case is realized for mutually orthogonal sets of 
basis functions fj n and /JJ where S mix = 0; if, in ad- 
dition, detT mix = 0, we have x mix (A) = (see also 



Sec. HHC] below), else x mix ( A) = r mix (e lA - 1)^ with 
r mix = det T mix /VdetS/'S/ n . 

Second, let us assume that S mlx is invertible, 
detS mlx 7^ 0. Let t™ 1x be the roots of the polynomial 



det[T" 



*] = 0. 



(33) 



The matrix T mlx (S mlx ) _1 then can be brought into a Jor- 
dan canonical form with t™ 1x on the diagonal and the 
characteristic function assumes the simple form 



Hpt ^ mlx „ 

dctb TT (1 

VdetS/'S^i 1 ' 



rr ix e iA ) 



(34) 



The procedure outlined above is straightforwardly gen- 
eralized to states with higher Slater rank. 



III. TWO PARTICLES 

A. Full counting statistics 

The above findings have interesting generic conse- 
quences for the charge transport of fermionic particles; 
in the following, we discuss the simplest case of two par- 
ticles, see Fig. [TJ where nontrivial exchange properties 
manifest themselves. For N = 2 particles the diagonal- 
ization (|24| can be carried out explicitly for arbitrary 
matrices T* and . The two eigenvalues n 2 are given 
by 



Tl,2 



=F %/a 2 - det If det Sf 
detS? : 



(35) 



where the parameter 2a = S^T^+Sf-T^ — 2Re(S{ 2 T2 1 ). 
Alternatively, the eigenvalues < r m < 1 are given by a 
minimum/maximum property^ 2 - 



ri = min T(g), 

geH 2 \S(g) = l 



T2 = Tti)> (36) 

geH 2 \s( g )=i 



with the eigenvectors 51,2 (&) given by those functions 
where the minimum/maximum values are attained, i.e, 
1X51,2) = Tl,2- Once the eigenvalues r m are known, the 
characteristic function X2 assumes the simple generalized 
binomial form 



X 2 (A) = (1 - n + ne lA )(l - r 2 + r 2 e lA ) 



(37) 



As a result, we find that in the new basis g m , the two 
particles traverse the scatterer independent of one an- 
other, i.e., the characteristic function is a simple prod- 
uct of independent one-particle characteristic functions. 
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Even more, the characteristic function is determined by 
the Hubert space spanned by the incoming states /i )2 and 
is independent of the choice of basis. Exchange effects 
manifest themselves when comparing the result (|37[) for 
the Slater determinant oc det/ m (fc„) with the result 
X2 ist (A) = (l-Tfx+T^Xl-T^+T^) for distin- 
guishable particles, \I> dlst cx /i(fci)/2(fe2): Exchange ef- 
fects are absent if both matrix elements S 21 = (/2I/1) = 
and T|i = 0, i.e., for orthogonal initial and transmitted 
states. On the other hand, a finite overlap of at least one 
pair of these states generates finite exchange effects via 
the substitution of T( nm in X2 lst by the eigenvalues r m in 

X2- 

The minimum/maximum property described above en- 
tails a set of a priori inequalities for the transmission 
probabilities P n involving the transmission matrix ele- 
ments T min = min{T{ 1 , J f 22 \ and T max = max{T{ 1 , T{ 2 }\ 
note that while the probabilities P„ do account for ex- 
change effects, the single particle matrix elements T( nm 
obviously do not. With initial (nonorthogonal) wave 
packets f m normalized to unity, S(f m ) — 1, the search 
for the extrema in Eq. (|36|) includes these states as 
well. We then obtain the set of inequalities < t\ < 
Tmin < 7m ax < T2 < 1. Using them to estimate 
P = (1 - n)(l - r 2 ) < (1 - r 2 ), P 2 = t x t 2 < n, and 
P x = 1 — P — P 2 , we can derive the following bounds 

Pa < 1 — ?max, Pi > Tmax ~ Pmin, P 2 < Tinin- (38) 

for the transmission probabilities for two particles. The 
above bounds set an upper limit on bunching (P 2 and Pa) 
and a lower limit on anti-bunching (Pi). Note though, 
that the bound on P 2 does not exclude an increase (due 
to exchange) of the transmission probability beyond the 
'classical' value P^ lst = T^T^ for distinguishable par- 
ticles, see X2 lst above. Indeed, since T^T^ < Pmin, a 
value P 2 ^> T{ 1 T22 remains possible. Such a result has 
been recently observed:^ the probability of two-electron 
events in the electron emission from a CsaSb photo- 
cathode in a photomultiplicr tube has been found to be 
much larger than the square of the probability for single- 
electron emission. This was observed both in the case 
of thermal emission without photocathode illumination 
and photoemission under weak photocathode illumina- 
tion. Furthermore, as detailed calculation shows, a large 
P 2 can also be obtained for wave packets with amplitudes 
f 2 (k) = fi(k + 5k) shifted in fc-space and a large over- 
lap integral S 21 , combined with a transmission amplitude 
suppressing k- values in the overlap region. 

B. Restrictions due to binomial statistics 

An arbitrary two-particle scattering process is fully 
characterized by the three parameters Po,Pi,P2, from 
which only two are independent; here, we assume that 
we can transmit only integer charges (no charge fraction- 
alization). In Figs.^a) and (b), we find the regions with 



different statistical properties that can be generated in 
a two-fermion scattering process, both in P$-P 2 param- 
eter space as well as in the noise {(n 2 )) versus average 
number (n) diagram. We start with the definition of the 
physically accessible regime in these diagrams: requiring 
that Pi = I-P0-P2 > (Fig.da)) andP , Pi, P 2 > 
(Fig. 0(b)), we find that the black regions are forbidden. 

Traditionally, starting from Poissonian statistics (P = 
1) relevant for the coherent light emitted from a laser or 
for the transport of a classical electron gas in a vacuum 
tube, much emphasis has been put on the distinction be- 
tween sub- and super-Poissonian statistics, with reduced 
and enhanced noise intensity as quantified by Fano fac- 
tors F < 1 (sub-Poissonian noise) and P > 1 (super- 
Poissonian processes). It appears to us that in the con- 
text of degenerate fermions, the generic starting point is 
the binomial statistics, instead, and more relevant qual- 
ifications are given by the regimes of sub-binomial and 
super-binomial processes introduced below. 

Nevertheless, let us start our analysis with the tradi- 
tional classification comparing a process with Poissonian 
statistics, which is realized on the dotted line in Fig. [2a) 
defined through the relation 

„ Po(l-Po) + P 2 (l-P2) + 2PoP 2 

F 1-P0+P2 =1 ' (39) 

i.e., 

P 2 = P - V2P0 - 1 with P > 1/2. (40) 

Within the dark-gray region noise is super-Poissonian, 
which is usually associated with the bunching of particles 
and therefore with bosonic statistics. Note that Fano- 
factors larger than the Poissonian value 1 require a large 
reflection probability Pq > 1/2; only when most of the 
particles are reflected one can observe the 'bunching' of 
the remaining transmitted objects. 

A much more natural classification for our fermion sys- 
tem is in terms of (deviations from) binomial statistics. 
The characteristic function % 2 for two fermions in a Slater 
determinant state can be cast into the generalized bino- 
mial form Eq. (|37p . which depends on two parameters 
Ti, r 2 . As a consequence, the probabilities satisfy the ad- 
ditional inequality 

V^+V^<1. (41) 

This condition follows from expressing the parameters 
Ti , t 2 through the probabilities Pq , P 2 using the relations 
Pq = (1 — 7~i)(l — r 2 ) and P 2 = tit 2 ; requiring a positive 
discriminant of the resulting (quadratic) equation implies 
the constraint (J41"jl which defines the light gray region in 
Fig. Wi a)), naturally termed the 'sub-binomial' regime. 
The (thick) black line bounding the general binomial (or 
sub-binomial) region is the line of usual binomial statis- 
tics, which is realized for the case of degenerate transmis- 
sion coefficients t\ — t 2 as they appear if the scattering 
does not depend on energy. 
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((n 2 )) 




unphysical 

super-Poissonian 
Poissonian (F = 1) 

super-binomial 
binomial 

generalized binomial 
(sub-binomial) 



FIG. 2: Diagrams describing the generic statistical properties of two-particle transmission, on the left as a P2-P0 probability 
diagram, on the right as a noise-charge {(n 2 ))-(n) diagram. The black regions are unphysical with probabilities Po,Pi,P2 
residing outside [0, 1]. The light gray regions describe generalized binomial (sub-binomial) processes, Eq. 11370. bounded by the 
black line characterizing usual binomial processes. The dotted lines correspond to a Fano factor F = ((n 2 ))/{n) equal to one. 
Within the dark gray regions noise is super-Poissonian with a Fano factor F > 1. Note that in order to observe super-Poissonian 
noise the reflection probability has to be large, such that Po > 1/2 and (n) < 1. 



The region with super-Poissonian noise (dark gray) 
and the sub- binomial region (light gray) are distinct, 
with the statistics of fermions incoming in a Slater- 
determinant state always residing in the sub-binomial do- 
main. Note that the counting statistics of an arbitrary 
two-particle process (without specification of exchange 
properties) also depends on two out of the three parame- 
ters Po, Pi, P2 (as the constraint Po + Pi + P2 = 1 needs 
to be fulfilled) but cannot be cast into the form Eq. (|3T|) 
in general, hence these processes are devoid of such an 
additional restriction. 

The P2-P0 diagram can be transcribed to the (experi- 
mentally more relevant) ((n 2 ))-(n) diagram, cf. Fig.[2^b). 
The physical constraints < Po, Pi, P2 lead to the set of 
inequalities, 

((n 2 ))>(n)(l-(n)), 
((n 2 ))<(n)(2-(n)), 

((n 2 ))>((n)-l)(2-(n)), (42) 

which can be cast into the more compact form (m + 
1 - (n))((n) - m) < ((n 2 )) < (n)(2 - (n)), with m = 
0, 1. The single large- and two small parabolas bounding 
the unphysical (black) regions are given by the second 
and the two (for m = 0, 1) first inequalities. For the 
generalized (or sub-) binomial statistics, the additional 
constraint assumes the form 

P=^<l-(n)/2, (43) 

with the equality applying to the binomial case with 
Ti = T2 . Within the gray region of the diagram the noise 
is sub-binomial F < 1 — (n)/2 and hence trivially sub- 
Poissonian, F < 1. Note that noiseless transmission of 



charge requires that an integer average charge is trans- 
mitted. 

The generalization of the above analysis to N incoming 
particles in a Slater determinant state is straightforward. 
The generalized binomial characteristic function is given 
by Eq. ©. The positivity of the probabilities P m > 0, 
m = 0, . . . , N imposes the N + 1 restrictions on the first 
two momenta (n) and ((n 2 )), (m + 1 — (n))((n) — m) < 
((n 2 )) < (n)(N ~ (n)), with m = 0, . . . , N - 1, defining 
a simple generalization of Fig. [DJb) with one large and 
./V small parabolas. In the generalized binomial case, the 
additional constraint 

F= -^-t^- <l-(n) /N<1 (44) 
(n) 

tells that the incoming Slater determinant states pro- 
duce a sub-binomial noise statistics. A similar result was 
found recently^! in the context of adiabatic pumping. 
The authors considered a time-dependent scattering ma- 
trix in the instant scattering approximation (i.e., an en- 
ergy independent scatterer) and obtained a generating 
function in a product form describing a generalized bino- 
mial statistics with parameters u m < 0; the u m relate to 
our r m via r m = (1 - w m ) _1 . 

C. Entangled states 

The above discussion for two particles lets us conclude 
that incoming Slater determinant states generate Fano 
factors F < 1 — (n)/2 < 1; such states are nonentangled. 
On the other hand, an entangled two-particle state can be 
generated with a sum of two Slater determinants; such 
an entangled state (with < a < 1, i.e., a state with 
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Slater rank 2, see Sec. Ill Fj) 

*(*i,Aa) = V^V I (k 1 ,k 2 ) + VT :r ^¥ I (k 1 ,k2), (45) 

is sufficient to generate all possible types of two-particle 
statistics: we choose the Slater-determinant wave func- 
tions and 'i' 11 (incoming from the left) such that they 
occupy different parts of momentum space, e.g., , J /I has 
only components below k c and 'J' 11 above. Furthermore, 
let the transmission be T\ = T/ C <fe c below k c and T 2 — 
Tk>k c above. For such a setup, all the overlap integrals 
vanish, e.g., J(dk 1 dk 2 /4:Tr 2 )^ n *(k 1 ,k 2 )^ 1 (k 1 ,k2) = 0, 
and we obtain (cf. Eq. (|2"8]1) 



X 2W=a(l-T 1 +T 1 eS x y + (l-a)(l-T 2 +T 2 e i y, (46) 

that is, the generating function is simply the weighted 
sum of the two individual generating functions for the 
Slater-determinant states. The statistics of such entan- 
gled wave functions is described by points in the P 2 -Pq 
diagram of Fig. [DJa) which lie on a straight line between 
the point p 1 for 'J' 1 and the point p u for 'J' 11 with a pa- 
rameterizing the line. Both p l and p u are situated on the 
binomial line, while the line connecting them may enter 
the super-binomial or even the super-Poissonian region: 
for example setting T\ — and T 2 = 1, the characteristic 
function is given by xi = a + (1 ~ a)e 2lX and F = 2a, 
which assumes values between zero and two (note that in 
the limit a — *■ 1, the wave function (|45[) is of Slater rank 
1, but nevertheless, the Fano factor approaches F = 2. 
As the Fano factor for Pq = 1 assumes the form 0/0 its 
value depends on the direction from which Pq = 1 is ap- 
proached). As simple Slater determinants produce only 
Fano factors up to 1 — (n)/2, a larger value serves as a 
test for the entanglement of the two particles . 35 ' 36 . For 
TV incoming particles in an entangled state of rank 2, the 
analogous construction (cf. Ill F[) produces a Fano factor 
F = Na with < a < 1, i.e., super-Poissonian statistics 
can be admitted for sufficiently large a. 



D. Two spin 1/2 particles 

Next, we consider the situation in the setup of Fig. [1] 
with incoming particles in normalized states fi(k) and 
f 2 (k) with overlap S = S 21 = (/2I/1) and carrying 
a spin 1/2 degree of freedom. We consider the case 
of spin-independent scattering, hence the coefficients in 
T* depend exclusively on fi(k) and f 2 (k). The four 
properly symmetrized states available to the two incom- 
ing particles are denoted by \]/ Sima (k), with s = the 
singlet {m s — 0) state and s = 1 the three (m s — 
— 1,0, +1) triplet states. The degrees of freedom k in- 
volve the momenta k m and spins s m of the particles, 
k = (ki, si; k 2 , s 2 ). The triplet states with m s = ±1 
are simple Slater determinant states 

*i,±i(k) = f= 1 ,^, n = [A( fc i)Xt/l( s i)/2(fc2)xt/l( s 2) 



- [(fcl,Si) <-> (fc 2 ,s 2 )]]. 



The characteristic function \i f° r the full counting statis- 
tics then is of the generalized binomial form with Tx/ 2 
given by Eq. (|35|) . 

Xi,±i(A) = (1 - n + r ie lA )(l - r 2 + r 2 e 4A ) (48) 
_ (1-T{ 1+ T{ ie '*)(l-T£ 2 +T£ 2 e'*) 

The states with m B = are more interesting as they 
are of Slater rank 2. Defining 

fl(k,s) = fi(k)xi(s), f 2 (k,s) = h(k)xi(s), 
/ 1 II (fc, S ) = / 1 (fc) Xi ( S ), /i I (fc, S ) = / 2 (fc)x T W, (49) 



we have 

*o/i,o(k) 



[f J (k) =F *"( k )] (50) 



with \f ,I / II (k) the normalized two-particle Slater deter- 
minants made from the states fin. The calculation of 
the characteristic function follows the procedure outlined 
above: As the matrices fiU = T^Smn and S™ = S mn 
are diagonal (the particles 1 and 2 are distinguishable), 
we immediately have 

X 1/U (X) = (1 - Tf a + T( ie * A )(l - T{ 2 + T{ 2 e iA ). (51) 

For the calculation of x mlx (A), the matrices T mix and 
S mix need to be evaluated. In the present case, they are 
purely off-diagonal with the off-diagonal matrix element 
given by T 2 y x = T 21 and S™* = S. Calculating the de- 
terminants in Eq. (|29p . we obtain the mixed component 
in the form 

X mix (A) = -(S~J f 21 +T f 21 e lX )(S*-J f 12 + J f 12 e^) (52) 
and the characteristic function is given by 



Xo/i,o(A) 



(l-T{ 1+ T{ ie « A )(l-T/ 2 +T 2 V A ) 



l±\S\< 



(S -n.+T^jS* ~1{ 2 +T[ 2 e^) 
1±\S\ 2 



± 



-f 



(53) 



(47) 



The result ([53]) agrees with the results in Ref. [l^. The 
characteristic functions for the two spin triplet states 
s = 1 with maximal magnetization m s — ±1 and the 
characteristic function for the triplet s = 1, m s = with 
zero magnetization coincide with the one for a Slater de- 
terminant of spinless fermions, Eq. (|3"T|) . This is because 
all three states involve identical orbital wave functions 
and the scattering process does not depend on the spin 
part of the wave function. The corresponding average 
number of particles (n)i ims and noise ((n 2 ))i ims reside 
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within the region of generalized binomial statistics, cf. 
Fig.! 

F hms < 1- <n)i,m./2. (54) 

The entangled singlet state (with s = 0) does not nec- 
essarily fulfill this condition. Rather opposite, for the 
case where the individual transmission probabilities of 
the two particles are equal, = T 22l the moments and 
the Fano factor always reside outside the region allowed 
by the generalized binomial statistics, 

F ,o > 1 - (n) ,o/2, (55) 

as a lengthy but straightforward calculation shows. 
Hence, this rather trivial setup can be used to discrimi- 
nate singlet from triplet states and also serves as an in- 
dicator of entanglement (as long as the inequalities are 
strict which is the case as long as S ^ and T 21 ^ ST^). 

A similar experiment was proposed by Burkard et al. 
Ref. fl8l which had two particles with equal energy come 
in from different arms in a symmetric beam splitter, see 
Ref. Qj] for a calculation of the full counting statistics for 
this setup. Our setup involves one single lead only, at the 
expense of requiring an energy dependent transmission 
probability (otherwise we end up on the binomial line 
which is devoid of any separation power). Furthermore, 
the discrimination between singlet and triplet states is 



determined by the presence or absence of generalized bi- 
nomial statistics and hence involves the binomial bound 
1 - (n) ,o/2 on F. 



IV. iV-PARTICLE TRAINS 

We consider the case of N incoming particles, all with 
the same shape of the wave function f(k) aligned regu- 
larly in real space with separation a; the wave function of 
the m-th particle then is given by f m (k) = f(k)e~ imak . 
The overlap and transmission matrices (|17[) are given by 
the Fourier transforms 

s f mn = J §U(fc)|V (m - )feQ , 

TL = J ^\f{k)\ 2 T k e i{m - n)ka ; (56) 

these are Toeplitz matrices as their elements depend only 
on the difference m — n between indices. In the limit 
N — > oo, the determinants of the Toeplitz matrices S* 
and — T* + T?e lX can be evaluated by reducing the 
integral over fc-space to an integral over the first Brillouin 
zone [0, 2n/a] and using Szego's theorem, see Ref. I20I and 

App.rjg 



logdetS'~iV / ^logj- V \f[(6 + 27Tm)/a}\ 2 }, 
logdet(S' -V + Ve iX ) ~ N [ ^log{- V \f[(0 + 2vm)/a]\ 2 [l - T [e+2mn)/a + T {g+27rm)/a eM. (57) 



The logarithm of these determinants scales linearly with 
N, a result that has to be expected as correlations be- 
tween particles vanish at large separation. Combining 
the results (1ST)) and replacing the integration over the 
angle 9 by an integration over the first Brillouin zone 
k G [0, 2n/a], we find the generating function in the form 

/" 27r / a rjh 

\og X N{\) = Na J Q ^ l0g(1 ~ n+ TkGlX) (58) 
with the effective scattering probabilities 

= Emgz l/( fc + 2irm/a)\ 2 T k+2lTm /a 

E mez l/(fc + W«)P ' 1 ; 

which denote transmission probabilities (with < r k < 
1) averaged over higher harmonics 27rm/a with weight 
|/(fc + 27rm/a)| 2 . 

Let us apply this result to wave packets generated by 
Lorentzian voltage pulses. As shown in Ref. [ToL a unit- 



flux (i.e., c J dtV(t) = hc/e = <!>o) Lorentzian voltage- 
pulse eVt (t) = 2fby/[(t — t ) 2 + 7 2 ], parametrized by its 
width 7 and time of appearance to, excites a single par- 
ticle with wave function f Xa (k) — ^47r^e~^ k ~ tXok €)(k) 
moving through the quantum wire (©(&) denotes the 
unit-step function; we remind that k is measured with 
respect to the Fermi momentum fc F ). Here e is the 
charge of the particle, xq — v F to parametrizes the po- 
sition, and £ = v F j the real-space width of the wave 
packet. A periodic sequence of unit-flux voltage-pulses 
V(t) = J2mez Vma/v F (t) applied to an interval to the left 
of the scatterer and driving one particle per time interval 
a/v F generates the transmission probabilities 

T k = (1 - e~ 4 </ a ) £ e- 4 ™«/ Q T fe+27rm/a . (60) 

m>0 

For nonoverlapping wave packets £ <C a, cf. Fig. [31a), 
exchange effects are absent. The sum in Eq. ([pTj]) be- 
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FIG. 3: (a) Train of nonoverlapping wave-packets, £ <C a. 
Each particle is transmitted independent of the others with 
a transmission probability f(dk/2iv)\f(k) \ 2 Tk depending only 
on its momentum distribution f(k). For wave-packets with 
width £, these probe transmission probabilities for momenta 
up to 1/^. (b) Train of strongly-overlapping wave-packets, 
£ 3> a. If the particles were transmitted independent of each 
other (no exchange effects), they would probe transmission 
probabilities Tk in the range up to l/£ ~ 0. Due to exchange 
effects, the particles fill up a Fermi sea determined by the 
density 1/a. Therefore, the particle train probes transmis- 
sion probabilities for momenta in the interval [0,27r/a]. This 
(stationary) state can be seen as a wave-packet analogue of 
the constant-voltage setup. 



comes an integral and the transmission probabilities as- 
sume the simple form r k = f(dk' /2n)\f(k')\ 2 T k , = (T) 
(as easily obtained from (|5§|) by replacing the sums with 
integrals), independent of k: Every particle probes the 
transmission probabilities with its weight |/(/c)| 2 . In the 
opposite limit, £ ^> a, i.e., for flat wave-packets which are 
strongly overlapping, the result for distinguishable par- 
ticles is r^ lst — > To. However, exchange effects force the 
system to fill the first Brillouin zone k S [0, 2%/ a] and 
the particles probe the transmission within this energy 
interval r k — T k (as obtained from (|60[) retaining only 
the m = term), cf. Fig. [21b). Taking the limit £ — > 00 
corresponds to the case of a constant applied voltage V 
of magnitude eV = 2irhv F /a (we remind that the time 
integral over one voltage pulse generates one flux unit 
$0 = hc/e and pulses are separated in time by a/v F ) and 
the generating function assumes the form 



log X ;v(A) = 7V- 



2nhv v r v/hvw dk 



eV 



2tt 



log(l-T k + T k e iX ). 



(61) 

This result then is the characteristic function for the full 
counting statistics for a constant voltage V applied to the 
left of the scatterer including an energy dependent scat- 
terer. The equivalent result has been found in Ref. [l^ if 
we perform the ad hoc replacement of the particle number 
N by a 'measuring time' t, N — > tv F /a — teV/2Trh (note, 



although this replacement appears sensible, it is non rig- 
orous as we have assumed the limit t — > 00 within the 
present scattering formalism; we will further comment 
on this later). 

The next order term in the asymptotic expansion 
for N — > 00 can be obtained using the generaliza- 
tion of Szego's theorem (Fisher-Hartwig conjecture, cf. 
Eq. (|A10P ) and is given by 



AlogXAr(A) = 



log TV 



4tt 2 



l-T, 



2-k/o, ' 



1-Tn 



(62) 



the logarithmic nature of this correction is due to the 
energy dependence of the transmission coefficient T k , in 
particular, its jump To 7^ T 27r / a across the first Brillouin 
zone k € [0, 2~k/<x\ (for T = T 27r / a the correction term 
is of order unity, see (|A7[) ). The correction for the noise 
term is given by 



A«n 2 )) = 



(^2ir/o — To) 
2^ 



■log TV 



(63) 



and similar corrections are obtained for the third- and 
higher-order cumulants. 



V. GENERALIZATIONS 

A. Unitary evolution and time-dependent counting 

We want to generalize the generating function xn a s 
given by Eq. (Tl4|) to account both for the specific time- 
evolution of the scattering state and for different counting 
procedures. Throughout this discussion, it is convenient 
to apply the Dirac notation and we rewrite the Slater 
determinant (11311 in the form 



1 



fNl 



Bgn(7r)|^ (1) ) <8> • • • <8> \^ {N) ); (64) 



Eq. ([64]) describes the initial iV-particle wave function 
at time t = composed of orthonormalized one-particle 
states \4>m) (here, tt denotes an element of the permu- 
tation group Sjsr). The choice of orthonormalized wave 
packets is only for convenience: as seen in Sec. Ill D[ a 
Slater determinant is invariant under general linear com- 
bination of states it is composed of; in particular, an 
orthonormalized basis can be chosen. 
Let 



U = exp 



4 I ' dt'Hit') 



(65) 



be the unitary evolution operator generated by the sin- 
gle particle Hamiltonian Ti(t)£L In the absence of inter- 
action, the evolution of the total system is governed by 
the product operator Tn(U), where, given a one-particle 
operator O, we define the iV-particle operator 



T N (0) = (D< 



O 



(66) 



JV-tij 
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acting simultaneously on all TV-particles. While we re- 
strict ourselves to noninteracting systems, we still allow 
for a time-dependent scattering potential which can gen- 
erate inelastic processes. The final state at time t is 
given by |^ ut) = rjv(W)|^). Including the counting 
field e ±zA / 2 , the wave- function assumes the form 



I*; 



IV(e ±iAS / 2 )|* c 



T N (e 



±i\Q/2 



LOW, (67) 



where Q is a projector (Q 2 — Q and — Q) on that 
part of the wave-function that has been counted. E.g., 
in the original setup of Ref. 4 with a spin at position xq 
and particles incoming from the left, the operator Qt = 
Jj dx \x}{x\ projects onto the causal interval I — [xq, Xo + 
v F t] (no such operator Q mimicking a spin-counter can 
be defined for particles incident from both sides); hence 
that part of the wave function which passed the counter 
during the time t picks up an additional phase e ±jA / 2 . 
Note that it is always the full phase A which is picked up, 
as the particle is either measured (eigenvalue 1 of Q) or 
not (eigenvalue of Q) . The characteristic function of the 
full counting statistics is given by the overlap (fidelity) 

Xn(X) = (*- t |*+ t ) - mr N (rfe aQ U)\y) (68) 

of the forward- and back-propagating wave-functions 
measured with opposite spin states. 

Next, we exploit that the expectation value of a prod- 
uct operator Tn(0) in a Slater-determinant state can be 
written as a determinant of one-particle matrix elements 
(4> m \0\(j) n } in the Hilbert space H N spanned by the states 



1 N 

mv N {pm = — j2 sgn(^ o To n (^ (m) |oi^ (m) ) 

7t,7t'G5jv m—1 

l N 

7T ,7T" €lSn m—1 

= det(0 m |O|0„); (69) 



this formula is at the origin of (most) results which cast 
the characteristic function of the full counting statistics 
into a determinant form. Making use of Eq. (|69|). we 
can rewrite the characteristic function Eq. (|6"8")1 as the 
determinant 

XN(\)=det(0 m \e lXutQU \0 n ) 

= det(<t> m \l-T Q +T Q e iX \<t> n ) (70) 



with 



Tn = U^QU; 



(71) 



in going from the first to the second line in Eq. (|70| . the 
exponential has been expanded and use has been made 
of the fact that Tq is a projector. With Tq a projector 
in the one particle Hilbert space H, its eigenvalues in the 
subspace lie between and 1 and Eq. (170)) leads to 
a generalized binomial statistics. 



In order to familiarize us with this new formula, we 
reproduce the results of the above section. We then are 
interested in the situation where the initial state is 
localized to the left of the scattering region and the final 
state describes the t — > oo asymptotic behavior where all 
particles have completed the scattering process. Within 
the basis of states left/right of the scatterer with momen- 
tum k 7 the asymptotic form of the propagator is given by 
the unitary (scattering) matrix 



UP. 



(72) 



where the coefficients (r' k ) and t k (t' k ) are the re- 
flection and transmission amplitudes of a particle in- 
coming from the left (right). The total propagator as- 
sumes the form U °° — § {dk/2Tt)\k) on t~U k X} - m {k\ where 
we have introduced the asymptotic states |fc); n ( out ) = 
(|fc)L,in(out), |fc) R ,i„(out)) which are in- (out-) going plane 
waves in the left/right lead: in a formal derivation, we 
have to consider the t — > oo limit of the evolution in (|6"5)) 
within an interaction picture with a trivial reference dy- 
namics U = J(dk/2Tr)e~^ kt (\k) inin (k\ + \k) out out (k\). 
The counting operator Q is given by the projection on 
the right outgoing lead, Qr, = (0, 1)^(0, 1), and we obtain 



r Q °° R = Hdk/2ir)\k) ta (t k y k )^t k y k )i 



(73) 



Since the initial single-particle wave functions i n {k\4> m ) = 
((k\4> m },0) are located to the left of the scatterer, the 
characteristic function assumes the form 



X (X) =dct(0 m |l-r + Te a |<M 



(74) 



with T = J(dk/2ir)T k \k){k\, in agreement with (Til)) ; 
here, we have shortened the notation \k) = \k)i n ^ in 
agreement with previous sections. The generalization of 
the result Eq. ([71)1 to many channels is straightforward: 
the propagator Eq. ([72"]) exhibits a block structure with 
matrices t k and r k describing the transmission and re- 
flection in the channel basis, the transmission probabili- 
ties Tfc = t k t k assume a matrix form and the state vec- 
tor \4> m ) adopts an additional channel index. Assuming 
an implicit summation over channel indices, the form of 
Eq. (|74p remains unchanged. The same comment holds 
for the spin index. 



B. Density matrix — finite temperatures 

The determinant in Eq. ([70")) is restricted to the sub- 
space spanned by the initial states \4> m ). Introducing 
the projection operator V = X)m=i l^roK^ml onto the 
subspace spanned by the initial states \4> m ), the deter- 
minant can be elevated to cover the whole Hilbert space. 
We split the total Hilbert space into the sector defined 
by the projector V and its complement projected onto 
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V± = 1 - V. The operator 1 - VT Q + VT Q e tX can be 
expressed in block form 



1)] 

l+T Q (e iA -l) T Q (eJ x -l) 




1 



(75) 



with the blocks operating in the VH and V±H sub- 
spaces. The determinant of the upper block-diagonal 
matrix Eq. (|75|) is given as the product of the determi- 
nant (JTUJ) in the P-block and the determinant of 1 in the 
Pi-block and thus the generating function assumes the 
form 



X n{\) = det(l - VT Q + VT Q e iX ), 



(76) 



where the determinant is taken over the entire one- 
particle Hilbert space H . 

Interestingly, this formula can be generalized to the 
case when the initial state is not a single Slater deter- 
minant, but an incoherent superposition of many Slater 
determinants with a density matrix of the form T(p)/Z 
in Fock space F = ® N H N , T(0) = ® N T N (0), Z = 
Tvp a T(p), and p is the one-particle density matrix, e.g., 
p = e~P( n ~^ for a thermal ensemble with temperature 
/? , chemical potential p and time- independent single- 
particlc Hamiltonian Ti ; here F a denotes the antisymmet- 
ric sector of the Fock space. Using the trace formula 38 , 
Trj? a [r(0)] = det(l + 0), where the determinant is over 
the one-particle Hilbert space, the characteristic function 
X(A) = Tr Fa [T(p)T(UU iXQ U)]/Z, cf. Eq. (JUD, assumes 
the form 



X (A) = det(l 

= det(l — r) + r\e 
= det(l — t)Tq + T]Tge iX ) 



pe lXU sw )/det(l + p) 



(77) 



with the one-particle occupation-number operator r\ = 
p/(l + p) (and arbitrary one-particle density matrix p); 
note again that the spectrum of tjTq resides between 
and 1 so that (|77|) denotes a generalized binomial statis- 
tics. 

As an example, consider the situation of two particles 
incident from the left, with wave functions 0i(fc) and 
02 (fc), where the process of particle generation is not de- 
terministic but involves some success probability: let pi 
(P2) be the probability that the first (second) particle 
is successfully created. In order to keep the discussion 
simple, we assume 0i(fc) and 02(&) to be orthonormal- 
ized (0 m |0n) = 5 mn . The initial state can be written 
as a density matrix Y(p)/Z with the one-particle density 
matrix 



Pi 



I -Pi 



P2 



1 ~P2 



102) (02 



(78) 



The weights in (|78|) are chosen to make the single particle 
occupation operator 77 = p/(l + p) have the form 



i.e., pi (P2) are the probabilities to occupy the state 1 
(2). The normalized density matrix [we use the fact that 
Z = Tr F T{p) =det(l + p) - 1/(1 - pi)(l - p 2 )] 

T(p)/Z = (1 - Pl )(l -pa) ®pi(l -p 2 )|0i)(0i| 

+ P 2 (1 -Pi)|0 2 )(0 2 | ©PlP 2 |0l)(0l| 8> |0 2 )(0 2 | (80) 

consists of three terms: The first term describes the 
zero particle sector which occurs with probability (1 — 
Pi)(l — P2). The second term involves one particle states: 
Pi(l— P2) [p2(l — Pi)] is the probability that only the first 
[second] particle is created. The third term shows that 
the probability to observe a Slater determinant of both 
states 4>i and 4>2 is P1P2; we have omitted terms which in- 
volve tensor products of more than one projector on the 
same state as they have no weight on the antisymmetric 
part of the Hilbert space. The generating function of full 
counting statistics is given by ([77| . 

X (A)=dct(l- ?? r Q 00 R +^e^) 

= (1 - piT + piTe lA )(l - P2 T + P2 Te lX ) (81) 

for the simplest case of asymptotic scattering with an 
energy-independent transmission probability, Tj, = T . 



VI. CONSTANT VOLTAGE 

Many results in the literature so far have been ob- 
tained in the stationary regime where a constant volt- 
age V is applied across the wire for long measuring times 
teV/h ^> i^ i 39 i 4 °i 41 Here, we discuss a wave- packet ana- 
log of the constant-voltage case. Contrary to the dis- 
cussion in Sec. IIVI involving a nonstationary finite train 
of N particles with the spin counter measuring all the 
time t — > 00, here, we consider a stationary situation 
in the thermodynamic limit (iV, L — > 00 with fixed den- 
sity n = N/L, L the system size) with two reservoirs 
disbalanced by the applied voltage V and the counting 
extending over a finite time t. 

We start with N particles residing in (left incident) 
scattering states 

<p k (x) = (e lfc * + r k e- ikx )&(-x) + t k e lkx Q(x), (82) 

with energies he = hv F k between E F and E F + eV . The 
scatterer is positioned at the origin. In order to regularize 
the problem, we go over to wave packets <j) m {x): We split 
the momentum interval [0, eV/v F ] into compartments of 
width Kk — eV/v F N and define the weights 



j, I y/2Tr/x, x(m — 1) < k < xm, 

Id, elsewhere, 



(83) 



with m € {1, . . . , N}. With the real weights f m (k), the 
(normalized) wave packets 



T] =Pl|0l)(0l| +P2|02)(02 



(79) 



0rn(a;) 



Y fm(k)fk(x) 



(84) 
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define states centered around the origin. Note that 
adding arbitrary global phases to the wave packets </> m (x) 
does not change their Slater determinant (up to a triv- 
ial global phase of the many-body wave function, see 
([64]) ). Keeping V constant and letting x — ► 0, the wave 
packets spread out in real space, the particle number 
N goes to infinity, the homogeneous particle density as- 
sumes the finite value eV/2nhv F , and the resulting cur- 
rent (J) = (e/h)TV is constant in time, cf. (fTOTj) . This 
procedure then properly emulates the constant voltage 
setup, as it generates the identical zero temperature den- 
sity matrix as the one obtained in a second quantization 
formulation by filling scattering states within the interval 
of width eV. 

In making use of the expression (|70p . we need the time 
evolution of the wave packets, 



d k 

<t>m(x;t)= I —e-^ kt f m (k)ip k (x), 



(85) 



as well as the counting operator Q t = Jj dx\x)(x\ pro- 
jecting particles on the space interval / = [xo,xo + v F t], 
where we assume the counter to be placed to the right of 
the origin, xq > 0. 



A. Generalized binomial statistics 

The characteristic function Xt(A) is the determinant of 
the matrix, cf. Eq. (f70ll . 



e autQ t u^ n) = ^ m{ t)\ e ^\cl, n (t)) (86) 
= J dx<t> m (x;t)*(x\e iXQ *\x)(j> n (x;t) 

<^mn ~t~ (c" 1 ) Qmn •> 



Xt (A)=dct[<5„ m + (e iA -l)Q ri 



(87) 



with 



Qr. 



dk'dk 
~4^~ 



t k f n (k)K t (k-k')tl,f^(k') (88) 



and the kernel 

pX + v F t piqxofl _ p -iqv F t\ 

K t {q)= dxe^ x - v ^ = - ^— ^ '- 

Jxo n 

= 2c i q ( Xo - VF tm sin(g« F t/2) 



(89) 



The matrix Q mn is a Hermitian matrix with real eigen- 
values < r m (t) < 1, hence the associated full counting 
statistics is generalized binomial for all times. The same 
result can be retrieved from Ref. [34] with an appropriate 
choice for the time-dependent scatterer. In the follow- 
ing, we discuss various limits for the generating function 
Xt(A). 



B. Short measuring time 

Assuming that N is large enough so that tk does not 
change appreciably over the interval x, i.e., xdktk 
1, the amplitude tk can be taken out of the integral in 
Eq. (|88|) . Assuming furthermore that the measurement 
time t is short, |<7|f F t < teV/h <C 1, we can expand K t (q) 
and obtain (to lowest order in teV/h) 



X f (A) = det S mn + (e lA - \)t*„ m t „ n xv F t 



x e 



3 (»_ ro) 2sin 2 (xa;o/2) 



✓2 ™2 



(90) 



The second term involves a matrix product 
(ui, 1)2, ■ ■ ■ , Vn)'(vi, V2, ■ ■ ■ , Vn) of a vector and its 



dual, where v m = t > 



*, and hence can be writ- 



ten as a projector, in Dirac notation, /i|w)(i;| with 
/j = 2(e lX — l)xv F t sin 2 (xxq / 2)/ it x 2 Xq. The determi- 
nant det(l + fi\v)(v\) then is given by the product of 
eigenvalues 1 + fJ-(v\v) (in the direction of \v)) and 1 (in 
the complement), det(l + /j,\v}{v\) = 1 + and we 

obtain 



X <<(\) = l + (e^ -l)xv F t 



2sin 2 (^a; /2) 

TTX 2 Xq 



N 



m— 1 



(«-0) 



l + a(e lX - 1), 



with (note that Nx = eV/hv F ) 



tv F 



eV/hVF dk, 
2^ 



T k < teV/2-Kh < 1 



(91) 



(92) 



(note that the dependence on the counter position xq has 
disappeared from the parameter a cx t in the constant 
voltage limit x — > 0). Pushing the calculation to higher 
order in teV/h, we find that in the expansion of Xt^(A) 
both, second- and third-order terms, vanish and the next 
correction appears only in fourth order, 

Ax f (A) = { - X ~ 1)2(<Wf)4 f eV/hVFdkdk ' ™ «. u» 
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-T k T k ,{k-k'Y 
(2tt) 2 



< 



te" - 1\ 2 (teVy 



(93) 



hence, for short measuring times the majority of counts 
involve either no or a single particle, while the obser- 
vation of two-particle events P2 < (teV / K) 4 / (2An) 2 is 
strongly suppressed, a consequence of the Pauli exclu- 
sion principle. Note that in the short measuring time 
limit, the specific nature of the counting device matters. 
Above, we have assumed that all intrinsic timescales of 
the counter are much shorter than the measuring time. 
Furthermore, we have neglected the effect of the Fermi 
sea which will produce an additional contribution. Nev- 
ertheless, even modeling the counter more realistically, 
the Pauli principle with its reduction of two- and more- 
particle events is expected to reveal itself. Furthermore, 
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it is possible to realize an experiment where the effect of 
the additional Fermi sea is absent: by applying a voltage 
to a quantum wire which is larger than the Fermi energy, 
particles incident from the right are blocked by the band 
bottom and only left going states within an energy inter- 
val E F (replacing the bias eV) contribute to the particle 
current i^ 2 . 

Note that the generalized binomial statistics Eq. (|3|) 
reduces to the simple Poissonian result 



logx(A) 



(e^ - 1) 



(94) 



in the limit of small generalized transmission probabil- 
ities r m <C 1; the result then only depends on one pa- 
rameter Y] m T m = tr7g t . In the limit of short measur- 
ing times, the smallness of the transmission eigenvalues 
is imposed by the small space interval in the projection 
Qt and J2 m Tm ~ a > see Eq- dHJ)- The same result is 
obtained in the long time limit, see ([9"5]) . provided the 
transmission probabilities T k themselves are small. 



C. Large measuring times 

In the asymptotic limit of t — ► oo, the kernel Kt(q) 
ensures energy/momentum conservation, rendering the 
problem diagonal in the momentum basis J*^- However, 
adopting the t — > oo asymptotic limit is incompatible 
with a regular derivation of a finite result. Here, we con- 
sider instead the case of large but finite measuring time 
t, while adopting the limit of infinite particle number 
N — > oo when letting the width hx = eV/v F N go to zero 
at constant voltage V. 

In the limit TV — > oo, the characteristic function x*(-M, 
which is the determinant of the matrix in Eq. (|86[) . cf. 
Eq. ([TO]), is given by 



Xt(A) = det(l - VTQ t + VTQ t e lX ) 



(95) 



where V = J^ V ^ hVF (dk/2Tr)\k)(k\ is the projector on 
the subspace of occupied states; the form ([95]) can be 
obtained from (|86p introducing the projector V to ex- 
tend the determinant over the whole Hilbert space, cf. 
Eqs. ([75]) and (|TrJ|) . and using the determinant identity 
det(l + AB) = det(l + BA) to shuffle t£ to the left 
of K t , which itself is the momentum representation of 
the projector Q t . The expression Eq. ([9"5]) then corre- 
sponds to Eq. ([75]) with the substitution Tq = TQ t . As 
Qt is a projector, Q\ = Q t , we can rewrite Eq. ([95]) as 
det[l+(e lA -l)Q t 7'T Q t ]. This determinant only needs to 
be calculated in the subspace QtH as the matrix is unity 
in the complement. In the subspace QtH, we use the or- 
thonormal real-space (rather than fc-space, see Eq. ([83]) ) 
basis 



9i{x) 



1 / \fe, e(l — 1) < x — so < el, 



0. 



elsewhere, 



(96) 



with e = tvF/L the width of a real space segment and 
I € {1, . . . , L}, The matrix elements of VT assume the 
form 



(gi\VT\g m ) = 



eV/hVF dk 
2^ 



T k (gi\k)(k\g„ 



(97) 



^l^dk 4sin 2 (efc/2) i(l _ m)ke 
2ir k ek 2 



i.e., they form a Toeplitz matrix. Applying Szego's theo- 
rem and taking the limit of large t and L with e = v F t/L 
fixed but small (e <C hv F /eV, hence 4sin 2 (efc/2)/e/c 2 w e) 
we obtain the generating function 



logXt(A) = tv. 



reV/hv F dk 

Jo 2^ 



\og{l-T k + T k e lX ), (98) 



cf, App.rxj 

Using the generalization of Szego's theorem (Fisher- 
Hartwig conjecture, see (|A7I) and (|A10|) ^ 4 — . it is possible 
to calculate the next order term. As the argument of the 
logarithm (cf. (|Al0p . note that x{9) is to be replaced by 
En& p (e+2,n)/e[l + (e iX -l)T^ +27Tn)/e ] with P k = 1, k £ 
[0, eV/hv F ] and P k = otherwise) exhibits discontinuities 
at k — and k — eV/hv F , the correction to the leading 
term is given by the two contributions originating from 
the jumps at k — 0,eV/hv F , 



Alogxt(A) 



log(t/to) 



4tt 2 



]Tlog 2 [l + ( e a -l)T fc ], (99) 



k=0,eV/hv F 



with to some small time cutoff; this result leads to log- 
arithmic corrections for the second-order and all higher 
cumulants. For the noise, the correction is given 



A((n 2 )) = 



± eVjhv F 



2tt 2 



log(*/*o). 



(100) 



Here, the logarithmic corrections in Eqs. |99|) and (|100p 
are due to fluctuations in the number of particles in a 
finite interval of length v F t. Therefore, fluctuations do 
not disappear for T — 1 as in Eqs. (|62|) and ([63]) where 
the number of particles is fixed and noise stems only 
from partitioning. In addition to the noise originating 
from the voltage bias, there is an equilibrium contribu- 
tion due to the Fermi sea at any finite measuring time, 
cf. Sec. IVI Bl For asymptotically large times, the first 
contribution grows logarithmically in timei 2 - 

The reason why Szego's theorem is applicable to the 
matrix Eq. ([9T]) is the presence of time translation invari- 
ance: matrix elements between states localized at two 
different places/times depend only on their space/time 
separation and hence they form a Toeplitz matrix. The 
same reasoning does not apply to the momentum ba- 
sis and that is why we could not apply Szego's the- 
orem directly to the matrix in Eq. ([56"]) . The result 
Eq. (j98p (and its generalization to finite temperatures, 
cf. Eq. (|113p ) has been found by Schonhammer— us- 
ing a double projection in his counting procedure: in- 
stead of relying on Szego's theorem, use has been made 
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of the relation logdet(l + M) = trlog(l + M), followed 
by an expansion of the logarithm. Evaluating the trace 
of each term, the phase factors (see Eq. (f89| ) appear- 
ing in the cyclic product of the kernel K t cancel mu- 
tually. In the long-time limit, the kernels become di- 
agonal (ensuring energy-conservation) with one of them 
contributing a factor of i, thus rendering the cumulant 
generating function logxt(A) linear in t. In alternative 
approaches, use has been made of a mapping onto the 
Riemann-Hilbert problem 2 ^ (this procedure enables the 
calculation of the leading term as well as the logarithmic 
correction Eq. (fM| ). or of time periodicity 2 ^ introduced 
in order to render log Xt extensive in t (this way, only the 
term linear in t is obtained). 



D. Fano factor for intermediate regime 

In order to understand the crossover between the short 
and long time behavior of the carrier distribution, we 
calculate the Fano factor F and present the result as a 
function of v t — teV/2irh (the incident particle number 
during time t) in Fig. [Ha) for several values of the trans- 
mission coefficient T (for a scatterer with energy inde- 
pendent transmission). For small times, the distribution 
is Poissonian and hence F{v t — *■ 0) — ► 1. The binomial 
distribution valid at large times provides the asymptotics 
F(v t — > oo) — > (1 — T). In order to find the crossover in 
between, we determine the matrix Q (see Eq. ([88]) ) . 



n (*--p) ,* . 

Wmn * t ycm z xn' 



i(n— m)>t{xo — v-pt) 



sin[(n — m)xv F t/2] 
7f(n — m) 



in terms of which the characteristic function assumes 
the simple form ([57)1 and hence logxt(A) = trlog[<5 mn + 
(e lA — l)Q m „] (again, we consider the limit x — > 
at fixed voltage V). The average transmitted charge 
(n) = -idxlogxt\\=o, 

(n)=trQ = tWp / — T fc , (101) 
Jo 27r 

grows linearly with the measuring time t; the above re- 
sult coincides with those obtained from the short and 
long time expressions and The noise ((n 2 )) = 

— dflogxt\x=o assumes the form 



((n 2 )) =trQ 



trQ 2 

eV/hVF dk'dk 



-TyTi 



sin 2 [(/c' 



(102) 

k)v F t/2] 



{k 1 - kf 



(in the limit x — > considered here, both momenta do 
not depend on the position xq of the counter, as the wave 
packets are infinitely spread). In order to keep the anal- 
ysis simple, we assume an energy independent transmis- 
sion probability, = T, over the interval [0, eV/Hv F ]. 
The average charge then is given by 



The Fano factor F 
form 



((n 2 ))/(n) can be cast into the 



F = 1 - Tf{v t ) 



with 



fiyt) 



1 d x {i-\ x \) s ^p^ 



For small times i>t <SC 1, 

f(yt) = vt 
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(104) 



(105) 



(106) 



while / approaches unity in the long time limit !/ t » 1 
(7 w 0.5772 is Euler's constant), 



fiyt) = 1 - 



l0g(27TI/ t ) + 1+ 7 



ir 2 v t 



(107) 



The corrections to the simple binomial result produce a 
logarithmic in time increase of the noise ((n 2 )); the re- 
sult (|107[) coincides with Eq. (| 1 00[) for the case of energy 
independent scattering probabilities = T. This loga- 
rithmic dependence in the noise is due to the fluctuations 
in the number of electrons in a finite segment of the wire.— 
Analogously, the third cumulant ((n 3 )) can be calcu- 
lated; the (numerical) results, shown in Fig. 0Jb), in- 
terpolate between the Poissonian value {(n 3 )) / (n) = 1 
for short times and the binomial result ((n 3 ))/(n) — 
T(l — T)(l — 2T) for long measuring times. 



E. Finite temperature 

We consider the case where particles are emitted from 
a lead at finite temperature into vacuum, i.e., we assume 
a single Fermi reservoir of particles (incident from the 
left) which are scattered with energy-dependent trans- 
mission probabilities (to the right). At finite temper- 
ature, scattering states are occupied according to the 
Fermi-Dirac occupation as described by the one-particle 
operator 77 = [e^ 7 ^ - ^ + 1] _1 . The characteristic func- 
tion XtPO is given by Eq. ([77]) . with Tq = TQ t and 
the interval I = [xq, xq + v F t] defining the projector Q t , 
cf. Sec. IVI CI Following essentially the calculation in 
Sec. IVI C[ i.e., calculating the determinant in the basis 
(|96|) and applying Szego's theorem, the result 



dk 

logXt(A) = tv F I ^log 



l + T k n L (k)(e 



(108) 



(n) = TteV/2nh = Tv t . 



(103) 



with n h (k) = (k\r)\k) = [ e /3( to F fc -M) + can be ob- 
tained for fixed but long measurement times t. For 
high temperatures at constant particle density \J3 — > 0, 
nh(k) w e - ^' ftvp ' !_ w] ) all transmission eigenvalues Tfc — 
T k /(e^ hVFk -t^ + l) w T k e-^ hVFk -^ approach zero. The 
logarithm in Eq. (|108|) can be expanded and the emis- 
sion statistics for electrons leaving a Fermi reservoir in 
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FIG. 4: Fano factor F — ((n 2 ))/(n) in (a) and third moment {{n 3 }} / (n) in (b) for constant voltage with energy independent 
transmission probabilities T = 0.25, 0.5, 0.75, 1, as a function of the incident particle number v t = teV/2irh. Note that the 
Fano factor approaches the binomial value F — 1 — T for ^ > 1 whereas for Vt <C 1 it is always close to 1 irrespective of the 
transmission probability T. The third cumulant interpolates as a function of vt between the Poissonian value ((n 3 ))/(n) — 1, 
cf. Eq. l[94"|l. and the binomial result ((n 3 ))/{n) = T(l - T)(l - 2T), cf. Eq. JSHJl. The oscillations (especially around i/ t fsl for 
T = 1) are due to the sharp edge in the occupation number at k — eV/hvF- 



the high temperature regime is given by a Poissonian 
statistics^ 

/rlk 
_ Tke -mvrk-») (109) 

The Fano factor assumes the value F = 1, independent 
of T fc . 

To complete the analysis, we discuss the extension of 
the constant-voltage result with two reservoirs to finite 
temperatures. We model the setup by two Fermi reser- 
voirs with occupation numbers Wl/rW = [exp[/3(fi,u F fc — 
Ml/r)] + f° r particles incoming from the left (L) or 
right (R) , respectively. The voltage enters via the bias of 
the chemical potentials eV = fit, — /ir. 

Unfortunately, it is not possible to define a projection 
operator Q which acts after the evolution and which can 
emulate the action of the spin counter, cf. the discussion 
below Eq. (|57|) . The reason is that there is no way to tell 
for a particle outgoing to the left at time t — > oo whether 
it was coming in from the left and was reflected at the 
scatterer (hence, no counting is done with the counter 
to the right of the scatterer) or whether it was coming 
in from the right and has been transmitted through the 
scatterer (hence passing the spin counter once). One so- 
lution to this problem is to perform a first projective 
measurement at the initial time j 27 i 29 this corresponds to 
replacing Eq. (|TT[) by the expression 



det(l + pWe tXQt Ue- tXQt ) 
det(l + p) ' 



(110) 



with the occupation-operator r\ = p/(l + p) = 
J(dk/2n)\k) in diag[n L (fc), n R (fc)] i n (fc|, the single-particle 
evolution U involving the scatterer but not the spin 
counter, cf. Eq. (|7^|) . and Q t a projector emulating the 



counting measurement of transmitted particles via pro- 
jection of the wave functions onto the lead to the right of 
the scatterer. The additional factor exp(— iXQt), as com- 
pared to l|77p , corresponds to the additional measurement 
before the evolution. 

In this article, we want to stick with the spin counter 
as a measurement apparatus. Contrary to the situation 
in Sec. |V] the action of the spin counter cannot be mod- 
eled by a projection onto the outgoing states, i.e., the 
operators U± do not separate any more into factors de- 
scribing the scatterer and the counter, U± ^ e ±lA ^*/ 2 £Y. 
Therefore, we have to make use of the full evolution op- 
erators U±, l^out) = r(W±)|*&), in the presence of both 
the scatterer and the spin counter, where the index '±' 
refers to the two spin states of the counter. The overall 
evolution (cf. Fig. [5]) then can be written as 



(111) 



where U, cf. Eq. (f72|) , is the evolution without account- 
ing for the presence of the spin-counter. The generating 
function for the full counting statistics assumes the form 

det(l + pe- lXQ ^ 2 U^e tXQ tUe- tXQ t/ 2 ) , . 

XtW = d^uTTp) ■ ( 1 } 



The two counting procedures Eqs. (| 1 1 0|) and (|112|) agree 
if the particles are only incident from the left, as the addi- 
tional counting factors, compared to Eq. (|77|) . contribute 
unity. For particles incoming from both left and right the 
two counting procedures do not necessarily coincide; only 
if Qt commutes with p, we can shift the factor e~ lA2 */ 2 to 
the left of p and then cyclically permute the factors in the 
second term of the determinant to assert the equivalence 
of (JIIDJI and (fm]). 
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-i\/2 t -> — OO 



VII. CONCLUSION 



JX/2 



out 



iA/2 



t — > +00 



FIG. 5: Sketch of incoming (t — » —00, above the dashed line) 
and scattered (large gray box) outgoing (t — » +00, below the 
dashed line) states measured by a spin counter (small gray 
boxes) placed to the right of the scatterer. Left-incoming and 
scattered states are described by black arrows, right- incoming 
and scattered states correspond to grey arrows. The action of 
the spin counter (in the up state) is included in the expression 
U+, see Eq. The evolution W+ involves two projections, 

one described by exp(— iXQt/2) at time t — > —00 and the 
second exp(-HAQt/2) at time t — » 00. The right (grey arrow) 
incoming state at t —* -co acquires a phase factor exp(— iA/2) 
(gray box on the right) in the first projection; this provides 
the correct counting field for its transmitted part (gray box on 
the left). The phase factor of the reflected part is canceled by 
the second counting operator (crossed box on the right). The 
left incoming state is unaffected by the first counting and its 
transmitted part acquires a phase exp(iA/2) at time t — > 00. 



The interpretation of Eq. (|112[) as the generating func- 
tion for the full counting statistics using the spin counting 
procedure faces problems, since Eq. (|111[) is not neces- 
sarily 2tt periodic and hence the counting may involve 
a noninteger number of particles. In certain situations, 
however, the spin counter nevertheless leads to sensible 
results. In particular, for asymptotically long measur- 
ing times t — > co, the counting projection operator Q t 
becomes basically diagonal in the energy or momentum 
basis. Commuting cxp(— iAQ t /2) with p and repeating 
the calculation in Sec. lVI Cl i.e., calculating the determi- 
nant in the basis (|96[) and applying Szego's theorem, the 
result 



/dk r 
— log[l + T k {n h (k)[l - n R (fc)](e lA -l) 

+ n R (fc)[l-n L (fc)](e- a -l)}l, (113) 



We have used the first-quantized wave packet formal- 
ism to calculate the generating function XnW of full 
counting statistics of fermionic particles in various phys- 
ical situations, such as N particles incident in Slater de- 
terminant states of rank 1 (nonentangled) , rank 2 (en- 
tangled), or incoherent superpositions of Slater determi- 
nants in Fock space with undetermined particle num- 
ber. Our formalism captures various features such as 
energy dependent scattering probabilities as well as time- 
dependent scattering and time-dependent counting. 

We have presented our results in determinantal form, 
with further simplifications explicitly unveiling a general- 
ized binomial statistics in various cases. Applications of 
our results include a classification of possible statistical 
behavior of two-particle scattering events and a partic- 
ularly simple singlet-triplet and entanglement detector. 
In the context of coherent transport of noninteracting 
(degenerate) fermions, the natural reference point in the 
discussion of statistical properties is the binomial dis- 
tribution; energy dependent scattering naturally shifts 
the noise into the sub-binomial (or generalized binomial) 
regime, whereas additional correlations through entan- 
glement can generate super-binomial noise statistics. 

Our results, calculated at zero temperature, remain 
valid for /3 _1 <C Hv F /£,, i.e., sufficiently narrow wave 
packets with a small width £ in real space. Furthermore, 
we have calculated the generating function for the con- 
stant voltage case in the long-time limit for any temper- 
ature. For short measuring times our results are valid in 
the temperature regime <C eV and we have found a 
strong suppression of P n >2 due to Pauli blocking. 

The central element underlying the appearance of 
a (sub-)binomial statistics in fermionic systems is the 
absence of interparticle interactions and entanglement. 
This result remains valid for a time-dependent scattering 
potential and finite temperature. We have analyzed the 
modification introduced by entanglement and have found 
that super-binomial statistics may be generated. The in- 
clusion of interaction, particularly within the scatterer 
where interacting particles become entangled, remains an 
interesting open problem. 
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is obtained; alternatively, the result (|113|) can be ob- 
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APPENDIX A: (STRONG) SZEGO THEOREM 

The (strong) Szego theorem applies to Toeplitz ma- 
trices and reduces the calculation of the asymptotic be- 
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havior of their determinants to a simple integration (plus 
summation) problem. We define a Tocplitz matrix start- 
ing from a complex-valued periodic function a(9) with 
a(9 + 2n) = a(9). In addition, we require that its wind- 
ing number with respect to the origin is equal to zero. 
We define the Fourier coefficients 



/ 7r a ( e > e 

Jo 2?T 



(Al) 



and the associated N x N Toeplitz matrix with elements 

[Ajy(a)]m,n = a m -n (A2) 

depending only on the difference between the indices m 
and n (banded matrix), m,n — 1, ...,N. The strong 
form of the Szego theore m 21 ' 45 states that 

logdetAAr(a) (A3) 

OO 

~ N [loga] + ^2 ^[loga]„[loga]_ n 

ri=l 

asymptotically for N — ► oo, with 



[log a\ 



dd 
2^ 



\og{a(0)}e- m6 



(A4) 



the Fourier coefficients of log[a(0)]. The first term in 
(|A3|) scaling with N 1 is the result of Szego's theorem, 
while its strong form applies once the sum in the second 
term converges — this correction then scales with N° . 

Given the Toeplitz matrix X* = S f + (e iX - l)T f , cf., 
Eq. ([55]) . we show how to find its determinant Eq. ([57]) 
in the asymptotic limit of large N . Specifying the matrix 
elements 

/dk 
— |/(fc)| 2 (l -T k + Tk e a )e^ m -^ ka , (A5) 

we find the original periodic function x(9) by calculating 
the Fourier series 



<&) = 



= -J2 \f[(0 + 2nm)/c 

n * — * 



(A6) 



X [1 — T(6+27rm)/a + T(9+2nm)/a e% 



Note that, while the original function x(k) — |/(fc)| 2 (l — 
T k + ?fee lA ) was defined on the real axis, the new expres- 
sion n{k) = x(9 = ah) is restricted to the first Brillouin 



zone k e [0, 2ir/a\. Fourier transforming the logarithm 
of x{ff) according to Eq. ([ A4]) . we obtain the asymptotic 
expression for the determinant 



logdetX 



N 



de 

2^ 



\og[x{9)] 



y^n[loga:] w [loga:]_ w + o(l), (A7) 



consisting of a main term oc N , a first correction stay- 
ing constant as N — * oo, and a remaining correction o(l) 
vanishing as N — > oo. The (logarithm of the) determi- 
nant in Eq. (f5T)) is derived by setting T = in (|A6|) . 
Finally, we obtain the (log of the) characteristic func- 
tion by simple subtraction (we replace the angle 9 on the 
unit circle [0, 2ir] by k = 9/a in the first Brillouin zone 
[0, 27r/a]), to leading order in N 

r 27r / a dk 

\og XN (X) = Na — log(l -r k + r k e lX ) (A8) 

Jo 2vr 

with the effective scattering probabilities 



Tk 



£ meZ |/(fc + 2wm/a)\ 2 T k+2 

7rm / a 

E ma l/(H2m/a)P ■ 



(A9) 



For a function x(9) which is continuous on the unit 
circle, i.e., x(2ir) — x(0), the sum in (|A7[) converges and 
the corrections to (|A8I) are constant when N — > oo (and 
similar for s(9) = £ meZ |/[(^ + 27rm)/a]| 2 /a in the cal- 



culation of logdetSj^). A more subtle situation appears 
in the situation where x(9) and/or s(9) are discontinuous 
across the Brillouin zone (Fisher-Hartwig conjecture)^ 
This situation is the usual case as the wave function f(k) 
is discontinuous at the Fermi level k = 0. Thus, the sum 
in (|A7|) is divergent and the next term in the expansion 
of (US) scales with log A (cf. also Eq. ((Ml)), 



AlogXAr(A) = 
log^O+V^O-)] - log 2 [ s (0+)/ S (0-)] 



4tt 2 



log N, (A10) 



and is followed by a constant term. Here, the number 
of particles N is fixed and noise is due to partitioning; 
hence, the logarithmic corrections (|A10[) have to be at- 
tributed to partitioning (and not to fluctuations in the 
number of particles as for the constant voltage result 
Eq. (|99|) ). This is also consistent with the vanishing of 
the correction (|A10[) for T = 1 where x{&) = e lX s{9) and 
x(Q+)/x(Q-) = s(0+)/s(CT). 
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